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ABSTRACT 

We present a model describing the evolution of Fanaroff-Riley type I and II radio AGN, and the 
transition between these classes. We quantify galaxy environments using a semi-analytic galaxy 
formation model, and apply our model to a volume-limited low redshift (0.03 ^ z ^ 0.1) sample 
of observed AGN to determine the distribution of jet powers and active lifetimes at the present epoch. 
Radio sources in massive galaxies are found to remain active for longer, spend less time in the quiescent 
phase, and inject more energy into their hosts than their less massive counterparts. The jet power 
is independent of the host stellar mass within uncertainties, consistent with maintenance-mode AGN 
feedback paradigm. The environments of these AGN are in or close to long-term heating-cooling 
balance. We also examine the properties of high- and low-excitation radio galaxy sub-populations. 
The HERGs are younger than LERGs by an order of magnitude, whilst their jet powers are greater 
by a factor of four. The Eddington-scaled accretion rates and jet production efficiencies of these 
populations are consistent with LERGs being powered by radiatively inefficient advection dominated 
accretion flows (ADAFs), while HERGs are fed by a radiatively efficient accretion mechanism. 
Subject headings: galaxies: active, jets, evolution—radio continuum: galaxies 


1. INTRODUCTION 

The growth of galaxies is closely related to 
that of the supermass i ve b l ack holes at their 
centres (iMagorrian et al.l 119981 : lHaring fe Rixl 12004 
iGiiltekin et all 120091) . The activity of this central en¬ 
gine is linked to the evolution of the galaxy through 
infalling matter accreted onto the black hole powering 
the jets of the active galactic nuclei (AGN). The active 
nucleus can have a profound impact on its surround¬ 
ings and the host galaxy evolution, through either ra¬ 
diative or kinetic feedback. The kinetic mode feedback 
is dominant at the present epoch, whilst radiati ve feed¬ 
back is prevalent in the early universe (z ~ 2 - 3: lFabiaiil 
120121 ). Kinetic feedback has been invoked to explain 
the lack of star f orma tion in the most massive galax¬ 
ies dCroton et al.l 120061) . and the su ppression of coolin g 
flows in the cores of massive clusters (iFabian et al.ll2003lh 
The energy input from AGN heats up or expels the sur¬ 
rounding intracluster gas (ICM), decreasing the rate of 
accretion onto the central black hole until it is shut off 
completely. Accretion can resume once the gas cools, 
or is replenishe d throu gh _an i ntera ction or merger with 
another galaxy (jBahcall et al.lU997lh The increase in ac¬ 
cretion rate may eventually lead to the retriggering of 
the radio jets, completing the AGN feedback cycle. Evi¬ 
dence of the sporadic and intermittent nature of this feed¬ 
back process is found in ob servations of rising radio bub¬ 
bles (|Churazov et al.ll2001fl and sources wi t h multiple ra¬ 
dio lo be pairs (e.g. iGiovannini et al.ll 19981 iVenturi et al.l 
120041 . Quantifying the typical jet powers of these radio 
sources and the duration of their active state is crucial 
to understanding the energetics of AGN feedback. 

Radio-loud AGN consist of twin jets emanating from 
the active nucleus in opposite directions. These jets 
interact with the surrounding medium, leading to the 
formation of a diffuse emission region. The radio lu- 
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minosity of this region, termed a cocoon, arises due 
to the synchrotron emission of shocked electrons (or 
positrons) injected from the jet. Observationally, these 
radio sources are separated into two classes based on 
the di stribution of radio surf a ce brightness in th eir co¬ 
coons dFanaroff fc Rilevlll974) . iFanaroff k, Rilevl type II 
sources have well-defined jet termination shocks located 
towards the ends of edge-brightened cocoons, whilst type 
I sources have their regions of highest surface brightness 
close to the core. This morphological dichotomy has been 
argued to arise due to either the influence of extern a l en¬ 
vironmental factors on the jet structure (e. g.lLaing et al.l 
11994 iBicknelll 1 199.4 IKaiser fc Alexander! 119971b or bv 
parameters associated w i th the jet productio n mech¬ 
anism (e.g. iMeierl (2001 : iGarofalo et all 12010 b Sev¬ 
eral radio sources have been observe d to exhibit a 
mixed morphology (IKaiser fc Bcstll2007l h suggesting the 
FR-I/II dichotomy is at least partly due to the en¬ 
vironment. Both FR-I and FR-II radi o sources are 
thoug ht to initially expand supersonically dMarecki et al.l 
120031b and it has thus been suggested that FR-II ra¬ 
dio source s may evolve i nto FR-Is when their jets are 
disru pted (iGoual-Krishna fc Wiital 119881 : IKaiser fc Best! 
I2007f l. There are numerous models in the literature de¬ 
scribing t he temporal evolution of the powerful FR-II 
class (e.g. IKaiser fc Alexandedll997l : lBlundell et a.l.lll999t 

IKaiser fc Cotter! 120021b and some models for the FR- 
I morphological class (e.g. iLuo fc Sadlerl l2010l h These 
models typically make a number of simplifications to re¬ 
main analytically tractable. External density profiles are 
usually approximated by a single power law with a con¬ 
stant exponent, inconsistent with X-ray observations o f 
the gas density in clusters (e.g. IVikhlinin et al.l 120061 1. 
The separate FR-II and FR-I models also only consider 
the two limiting cases, where either the ram (for FR- 
IIs) or external pressure (FR-Is) component of the total 
pressure dominates. 

Properties of AGN central engines can be studied 
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using optical emission lines (see e.g. lHardcastle et al.1 
120091 : IBest fc Heckmanl I2012H . In this approach, AGN 
are partitioned into two populations, high- and low- 
excitation radio galaxies (or HERGs and LERGs), based 
on either line ratios or line equivalent widths. These 
emission line properties are closely correlated with ra¬ 
dio source morphological classification, with almost all 
FR-Is found in LER Gs and many, but not all, FR- 
IIs inhabiting HERGs (lHardcastle et al.ll2007t iLin et all 
l2?rmlh These two populations are thought to host black 
hole s powered by differen t accretion flow mecha nisms 
(e.g. lHardcastle et alJ2007t IBest fc Heckmanll2012l ). with 
LERGs f uelled by an advectio n dominated accretion flow 
(ADAF, iNaravan fc Yil 119951) and the radiatively effi- 
cicnt HERGs by a thin (iShakura fc Sunvaevl fl973h or 
possibly slim (jAbramowicz et al.lll988fl disk flow. The 
accretion flow state near a supermassive black hole is 
largely dependent on its accretion rate, though spin also 
contributes. Black holes with low Eddington-scaled ac¬ 
cretion rates m -C 0.01 are thought to be fuelled by 
ADAFs and those with higher rates of accretion b etween 
m = 0.01 and 0.3 powered by a thin disk (e.g. iMeierl 
12001b iPark fc OstrikeilfeOOll) . 

In this paper, we develop a hybrid model for both FR- 
II and FR-I type sources, and the transition between 
these morphologies (Sections [2] and [3j. We model the 
evolution of radio sources in gaseous atmospheres de¬ 
rived from a semi-analytic galaxy formation model and 
apply our model to a local AGN sample to characterise 
their energetics (Section B}. In Section 0 we examine 
the relationship between the radio cocoon-ambient gas 
interaction and radio source morphology. Finally, in Sec¬ 
tions [6] and 0 these results are compared to the theoreti¬ 
cal expectations for AGN feedback and used to examine 
the accretion properties of high- and low-excitation radio 
galaxies. 

The ACDM concordance cosmology with Qm = 0.3, 
Ha = 0.7 and Hq = 70kms _1 Mpc _1 is assumed 
throughout the paper. 

2. RADIO SOURCE MODEL 

In this section, we develop the dynamical model. 
Our model is large ly based on the FR-II model of 
iKaiser fc Alexanderl (119971). and the pre ssure-limited ex¬ 
pansion model of ~ lLuo fc Sadler |201(j ). We combine 
these two limiting cases in a single framework, and in¬ 
clude a more complete treatment of the environments 
into which radio sources expand. The key parameters 
defined in this section are summarised in Table 0 

2.1. Radio Source Dynamics 

The FR-II dynamical model comprises a relativistic 
plasma jet emanating ballistically from the active nu¬ 
cleus with constant opening angle (Figure [lji). A recon¬ 
finement shock will form further down the jet enabling it 
to transition from the ballistic regime into one of pressure 
equilibrium with its cocoon. The radius of the jet will 
remain constant after this reconfinement shock since the 
cocoo n pressure remains u niform (due to the high sound 
speed: iKaiser fc Alexander]1199 71. The jet terminates at 
an outwardly moving hotspot in a jet shock with the ram 
pressure of the jet material distributed over the working 
surface, generating a bow shock. The pressure exerted 
on the working surface is balanced by the pressure of 


the shocked ambient gas which surrounds the hotspot. 
The hotspot region is overpressured with respect to the 
rest of the cocoon, and backflow of overpressured plasma 
inflates the cocoon. The cocoon and bow s hock are ex¬ 
pected to expand ii i a se lf-similar manner ( Falld 119911 : 
IKaiser fc Alexanderl 119971) . The luminosity of the co¬ 
coon comes from the radio emission due to synchrotron 
electrons processing about the magnetic field lines of the 
cocoon. 

In the FR-I model (Figure [TJd), a relativistic jet 
emerges from the central black hole, inflating a cocoon, 
similar to the FR-II model. Particles will still be accel¬ 
erated by jet shocks, however, these shocks occur closer 
to the nucleus for the low-power jets of FR-Is. The re¬ 
gion where shocks occur is referred to as the flare re- 
gion, the FR-I cou nterpart of the FR-IIs hotspots (e.g. 
iLuo fc Sadlerl [2010 1. The diffuse emission region, which 
extends well beyond the shock, is regarded as the lobe or 
cocoon in analogy with FR-IIs. 

2.2. Environment assumptions 

We model the expansion of the pressure-inflated radio 
cocoon in a spherically symmetric external atmosphere 
with an arbitrary density and temperature profile. The 
external atmosphere is partitioned into small radial re¬ 
gions in which the external density and temperature are 
each approximated by a simple power law profile. The 
density profile in each small region is thus approximated 
by a power law of the form 

Px = Po(^) = fcr _/3 , (1) 

where r is the radial distance from the active nucleus at 
the centre of the source and the exponent /3 is constant 
in the region. The density parameter k is defined as 
k = po^cT for some radius ao, with corresponding density 
po, located in the region approximated by the relevant 
power law. 

The temperature of the external environment in each 
small region is similarly modelled as 


where the exponent £ is constant in the region, fcs is the 
Boltzmann constant, fh ~ 0.6?n p is the average particle 
number density for proton mass m p , and l is a constant 
of proportionality termed the temperature parameteiQ. 

2.3. Transonic dynamical model 

Dynamical models of powerful type II radio sources 
show that the overpressured cocoons expand self- 
similarly . i.e. t he shape of t he co coon remains fixed (Falld 
119911 : IKaiser fc Alexander! 1 199711 . As the cocoon expan¬ 
sion slows, however, the surface velocity of the inner 
parts of the cocoon will fall below the local sound speed 
of the ambient gas. The expansion rate in this subsonic 
phase (r oc tVG-ftb lLuo fc Sadleill2010h is lower than for 
supersonic expansion (r oc for realistic density 

1 The m/k b constant is introduced in this definition to simplify 
subsequent equations. 
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Figure 1. FR-II (left) and FR-I (right) model schematics. Jets and cocoons on only one side of the AGN are shown. 


profiles with /3 < 2: lKaiscr fc Alexander! 119971 ). and thus 
the assumption of self-similar growth cannot be extended 
to this transonic regime. We therefore model the cocoon 
as an ensemble of small volume elements in pressure equi¬ 
librium, where the growth of each element is otherwise 
independent. In the subsonic phase we further consider 
Rayleigh-Taylor instabilities expanding about the surface 
of this volume (see Section 0. 

In order to quantify the energetics of the adiabatic ex¬ 
pansion of the radio source we need expressions for the 
volume of the cocoon and its pressure. The evolution 
of the cocoon on one side of the central active nucleus 
is modelled explicitly, with the properties of the entire 
radio source calculated by assuming two identical, anti¬ 
parallel cocoons. Below we derive expressions for the size 
and luminosity evolution of the radio cocoons. 

2.3.1. Cocoon geometry 

The radio cocoon of our dynamical model is an en¬ 
semble of small angular volume elements in pressure 
equilibrium. Each element of fixed angular width is as¬ 
sumed to receive a constant fraction of the jet power as 
the cocoon expands. This assumption yields self-similar 
expansion when the whole cocoon is in the supersonic 
phase. The volume of each small angular cocoon ele¬ 
ment [9 — d6/2,9 + dd/2) is thus related to its radial 
length R(9) by 

dV ( 9 ) = ^ sin 9d9, (3) 

o 

where 9 is the angle between the jet axis and the radial 
line passing through the volume element (Figure [5]). The 
total volume of the cocoon on one side of the central 
nucleus is the volume integral of Equation [3] over the 
domain 9 € (0,|). To provide initial conditions for this 
model, we assume our cocoon is initially ellipsoidal in 
shape with a circular cross-section along the jet axis. The 
ratio of the length of the cocoon on one side of the central 
active nucleus, R(9 = 0), to the radius of the cocoon 
transverse to the jet axis is defined by A, the cocoon 
axis raticB 

2.3.2. Working surface pressure 

The pressure at the working surface of highly su¬ 
personic FR-II radio sources is normally calculated us¬ 
ing ram pressure arguments or the Rankine-Hugoniot 

2 Hence for a sphere we have A = 1. This axis ratio is rela ted 
to the aspect ratio Rt defined by IKaiser fc Alexander! (i 1997 1 as 
A = 2Rt- 


jump conditions for a plane parallel sho ck (e.g. 
IKaiser fc Alexanderl 119971 ILuo fc SadleH l2f)T(il b For a 
subsonic FR-I source, pressure due to such expansion 
is assumed insignificant with the working surface pres¬ 
sure taken to be equal to the external pressure (e.g. 
ILuo fc Sadlerll2010lh The working surface of the cocoon 
is defined here as the location where the expanding co¬ 
coon and external environment meet. 

For a compressible fluid and ignoring external forces, 
the working surface pressure is related to the external 
pressure p x {9) using the Rankine-Hugoniot jump condi¬ 
tions for a plane parallel shock. The working surface 
pressure approaches the external pressure as the cocoon 
expansion speed nears the sound speed, i.e. Mb —i 1 
where Mb ((9) is the Mach number of the bow shock. In 
the subsonic regime we do not expect the working surface 
pressure to fall below the external pressure. The pressure 
at the working surface should thus equal the external 
pressure for subsonic expansion; th is is the “pressure- 
limited expansion” of lLuo fc Sadlerl ( 2010 1. The general 
expression for the working surface pressure p{9) of each 
ensemble element in the cocoon is therefore 

2FxM r x +i rx ~ 1) Px for M b > 1 ^ 

Px for M b < 1, 

where T x = | is the adiabatic index of the exter¬ 
nal medium. This equation can be rewritten in terms 
of the radial length of the ensemble element R(9 ), and 
the density and temperature parameters, k and l respec¬ 
tively. The pressure of the external medium is given by 
Px = PxTxik-g,/fh) = (kl)R~^ + ^\ where the right hand 
equality is obtained using Equations Q] and [2] The Mach 



Figure 2. A cross section of the entire cocoon (solid line), with the 
non-synchrotron emitting, Rayleigh-Taylor mixed region shaded in 
light grey. Geometric parameters are shown in this schematic for 
the ensemble element at angle 0. The shaded rectangle of length 
h(t, 0) represents the height of the mixing layer. The height of the 
mixing layer increases towards the cocoon’s semi-minor axis where 
the normal velocity and thus Mach number are lower. 
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Table 1 

Environment and cocoon geometry parameters. 


Pressure profile parameters 


density parameter 
density exponent 

k 

9 

p x = kr & 

temperature parameter 

i 

t x = lr~t 

temperature exponent 

5 


Ensemble element parameters 


angle of element from jet axis 6 

fraction of jet power in element dX 

volume of ensemble element dV 

radius of ensemble element R 

radial velocity at element surface v 


Dimensionless constants 
dimensionless radius 77 = R{0)/R(6 = 0) 

dimensionless velocity £ = v±(6)/R(0 = 0 ) 


number of the bow shock is related to the expansion rate 
of the ensemble element as Mb = vj _/c x , where v±(9) = 
(R/p is the velocity normal to the surface (see Appendix 
for derivation and definitions of the dim ensionless radi us 
and velocity, rj and C) and c x (9) = -^Tx/cbt x (0)/to is 
the sound speed. The Mach number can then be writ¬ 
ten as Mb = [Rt((R/r]) 2 /(T^l)] 1 / 2 , and hence Equation 
0] becomes 

_ f ^kR-PiCR/r]) 2 - Z*=±(Jel)R-V + *'> for M b ^ 1 
P ~ \ (kl)R-V+® * for M b < 1 

(5) 

The average pressure inside the cocoon is calculated by 
weighting the surface pressure of each ensemble element 
by its volume, since the high sound speed is expected to 
yield uniform internal pressure. 

2.4. Adiabatic expansion 

The adiabatic expansion of each angular volume ele¬ 
ment in the cocoon ensemble is related to the pressure 
imparted on that element at the surface, p(9) 1 its volume 
dV ( 9 ) and the jet power Q by 

pdV + T c pdV= (T c -l)QdX, (6) 

where T c = | (|Kaiser et a.1.1 Il997h is the adiabatic in¬ 
dex of the cocoon and the derivatives are with respect 
to the time t. The differential quantity d\(9) is the frac¬ 
tion of the jet power injected into the ensemble element 
[9 — d9/ 2, 9 + d9/ 2), derived in the Appendix. Note that 
Q is the kinetic power of one jet, i.e. the total AGN 
kinetic power is Q tot = 2 Q. The source linear size (an 
observable quantity) is then related to the length of the 
ensemble element along the jet axis by D = 2 R(9 = 0 ). 
The subsonic expansion phase, however, must be con¬ 
sidered separately from the supersonic phase due to the 
discontinuity in cocoon pressure in Equation [5] We note 
that the hotspot diameters of typical FR-II sources are 
approximately 100 ti mes smaller tha n the linear size of 
their cocoons (lHardcastle et al.l[l998T) . The energy stored 
in these high pressure hotspots can therefore be neglected 
from Equation [6] since it is small compared to that stored 


in the rest of the cocoon. 

2.4.1. Supersonic expansion 

We first consider the supersonic phase expansion of en¬ 
semble elements in the cocoon. Substituting the relevant 
expressions for the pressure and volume of each element 
into Equation [6] yields a second order differential equa¬ 
tion, which can not be solved analytically in general. We 
therefore adopt a numerical scheme using a fourth order 
Runge-Kutta method by rewriting the differential equa¬ 
tion as a system of two first order ODEs. Thus for each 
ensemble element in the cocoon we must solve the fol¬ 
lowing system of equations: 

R = v 

. = 3(r x + l)(r c -l)QR^- 3 dA (P - 3T c )v 2 

8ttv(C/ ri) 2 k sin 9d9 2 R (7) 

(T x — 1)(3T C — P — i)l 
4R«+!(C/r?) 2 

where v(9) is the radial velocity at which the cocoon 
expands along each angle 9. The first of these equations 
thus gives the velocity of the expanding cocoon and the 
second its acceleration. 

2.4.2. Subsonic expansion 

We now examine the expansion of ensemble elements 
in the subsonic phase. Substituting the expressions for 
the subsonic pressure (Equation [5j and volume (Equa¬ 
tion [3]) into Equation [5] yields a first order differential 
equation, which can be solved analytically. The solution 
to this differential equation is of the form R oc t 1 ^. 
However, this expression can only be used to model the 
evolution of a radio source that is subsonic for its en¬ 
tire lifetime. In particular, the ensemble element radius 
predicted by the subsonic solution at the sonic boundary 
differs from that of the supersonic phase (Equation [TJ). 
because the subsonic solution ignores the previous (su¬ 
personic) expansion history of the radio cocoon. The ex¬ 
pansion in the subsonic phase must therefore be solved 
numerically, using the supersonic phase solution as an 
initial condition. We thus rewrite our analytic solution 
as a system of differential equations, 

R = v 

. (P + £- 2)v 2 ( 8 ) 

V ~ R ’ 

where again the first equation gives the velocity and 
the second the acceleration for each ensemble element. 

2.4.3. Conditions of continuity 

The differential equations derived in the previous sec¬ 
tions are only applicable in an environment with constant 
density and temperature exponents (i.e. constant P and 
£). For general environments, in practice we approxi¬ 
mate the density and temperature profiles as a series of 
power law segments (typically 100 segments per source), 
and ensure the continuity of ensemble element sizes and 
pressures across the region boundaries. These require¬ 
ments are met by using the radius R and velocity R in 
each region as initial conditions for the next one. 
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Figure 3. Luminosity-size profile (or L-D track; red curve) 
for a 10 11 Mq stellar mass AGN host galaxy with a jet power of 
Qtot = 10 35,5 W in a typical environment. The transonic phase 
occurs when some parts of the cocoon are expanding subsonically, 
while others are still supersonic. The fall in luminosity beyond 
10 Myr is due to the Rayleigh-Taylor mixing of the cocoon with its 
environment. The dotted blue line shows the evolution predicted 
by our model using the supersonic solution throughout the source 
lifetime. Note the different time axes for the two models (lower axis 
for our model, and upper axis for the FR-II only model). These 
models yield very different sizes and luminosities at ages exceeding 
1 Myr. 


2.5. Radio source luminosity 

In order to observationally constrain the properties of 
radio AGN, we need to model the radio luminosity due to 
the synchrotron emission from the relativistic electrons 
in AGN cocoons. Synchrotron-emitting electrons will 
lose energy through adiabatic, synchrotro n and inverse 
Comp ton losses. We follow the formalism of lKaiser et al.1 
(|1997t ). modifying their Equation 16 to include our multi¬ 
power law pressure profile. The total cocoon luminosity 
(i.e. both lobes) at frequency v is then given by 




^ 42(1 r c )/r c 


Ve(t) 

V e (ti) 


-(r B +i/3) 


(9) 


dti, 


where or is the Thompson cross-section, c is the speed 
of light, Tb = 5/3 is the adiabatic index of the magnetic 
field, s = 2.14 is t he electron energy i njection index and 
9b p = (s + l)/4 dKaiser et al.lll997lh Here, 14 is the 
volume of an electron ensemble where 7 is the Lorentz 
factor of an electron emitting radiation predominantly at 
frequency i/, whilst 71 are the Lorentz factors at the times 
t\ when the electrons are initially injected into the radio 
cocoon. This Lorentz factor 7 ; is calculat ed iterative ly 
(over each power law) using Equation 10 of lKaiser et al.l 
iHH), and finally used to determine the number density 
coefficient no defined in their Equation 8 . 

We note that there are significant uncertainties in the 
estimation of the radio source luminosity at a given jet 
power since the physics o f these objects is not fully 
understood. IWillott et al.l (Il999f) described these un¬ 
certainties using a single factor /, which is observa¬ 
tionally constrained to lie between ~ 1 and 20. In 
Section [5j we show that our model is consistent with 


j et powers estimated fr om X-ray cavity measurements 
(|Heckman fe Best! I2014T) for / ~ 5. This factor flows 
through directly from our jet powers into our measures 
of the AGN energetics, but also subtly affects our source 
age estimates. For an uncertainty factor of / = 5, 
the estimated source age at a given size and luminos¬ 
ity is approximately 0.3 dex lower than that calculated 
using / = 1. That is, the source age decreases as 
roughly f -1 / 2 in ag r eemen t with Equations 4 and 5 of 
iKaiser fc Alexander (|1997f) for supersonic expansion. In 
this work, we calculate our results using / = 1 but retain 
this uncertainty factor in equations so these results can 
be readily scaled for any realistic choice of this variable, 
e.g. / < 5. 

A typical luminosity-size track produced by our model 
is shown in Figure [3j Our new model predicts a very 
different temporal evolution of AGN size and luminos¬ 
ity compared to the standard FR-II models. This has a 
large effect on the derived jet powers and ages of ob¬ 
served AGN. The luminosity increases initially in the 
core regions where the density of the host environment is 
roughly constant (with radius) before decreasing as the 
density profile steepens towards /? = 2. This behaviour 
is similar to that found in 3D M HD simulations of ra¬ 
dio je ts in a range of environments (lHardcastle fc Krausd 
120m By contrast, dynamical mod els which assume a 
single 0 ~ 2 power law profile (e.g. IKaiser et al.l 119971 : 
iBlundell et al.lll999l not shown in figure) predict the lu¬ 
minosity decreases monotonically with increasing size. 
Our model also predicts the source luminosity will rise 
suddenly upon entering the transonic expansion phase. 
This increase occurs since the cocoon pressure equili¬ 
brates with that of the host environment rather than 
continuin g to drop o ff at the r ate predicted by the FR- 
II model. iLuo fe Sadlerl (|2010f ) find a similar increase in 
luminosity using their pressure-limited expansion model 
for FR-Is. The fall in luminosity beyond 10 Myr is due to 
the Rayleigh-Taylor mixing of the cocoon with the host 
environment as detailed in Section [3] 


3. RAYLEIGH-TAYLOR MIXING 

As the radio source expansion slows down, the cocoon 
becomes susceptible to the Rayleigh-Taylor instability. 
Different parts of the cocoon become subsonic with re¬ 
spect to the ambient medium at different times. As the 
cocoon plasma becomes entrained by this surrounding 
medium, the high energy synchrotron electrons in the 
cocoon will be collisionally reduced to the ambient en¬ 
ergy level of the much denser surrounding atmosphere. 
The radio emission from these mixing regions is thus sig¬ 
nificantly reduced; here we assume it to be zero. The 
fluid mechanics of this situation are considered here in 
order to estimate the volume of the cocoon that remains 
unmixed, and thus visible, as the cocoon ages. Radio 
luminosities estimated using the equations of the previ¬ 
ous sections will be modified by this visible fraction to 
account for the Rayleigh-Taylor mixing. 

The expanding surface of the u nderdense cocoon will 
becom e Rayleigh-Taylor unstable (|Ravleighlfl8 83: Taylor 
Il95(ill as it enters the subsonic phase. The thick¬ 
ness of the Rayleigh-Taylor mixing layer can be calcu¬ 
lated based on the late time self-similar growth phase 
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(jFermi fc von Ncumannl 11953IFI . T he rate of e xpan¬ 
sion of this layer can be shown (iCook et al.1 12004 
iRistorcelli fc Clarkil2004D to be described by 


Here h is the thickness of the mixing layer (see Figure 
El), p x and p C oc are the densities of the external medium 
and cocoon respectively, and re is a dimensionless growth 
parameter. Most lab-based experiments have returned 
value s for re in the range 0.03 < re < 0.07 (e.g. iDimontd 
l2004f) . The effective gravitational acceleration, g e s, is 
the sum of the gravitational potential of the galaxy, 
Sgrav; and the acceleration of the cocoon, g coc = —v 
(with v as defined in Equation [8] for the subsonic phase). 
The gravitational acceleration on the cocoon is given by 
flgrav = c x 2 (/3/r x i?), where c x is the sound speed of the 
ambient medium with adiabatic index r x , and /3 is the 
exponent of the power law density profile j ust outside the 
cocoon (e.g. IShabala fc Alexandcrll2009al) . 

A general expression for the thickness of the mixing 
layer, for constant re but with a time-dependent g e s, is 
then given by 

h{t) = ^k 1/2 g eS 1/2 (t)dt + h 0 1/2 ^j . ( 11 ) 

The initial conditions to and ho ~ 0 in this equation 
are the time at which the flow first enters the self-similar 
growth phase and the corresponding thickness of the mix¬ 
ing layer. The integral in this equation can be solved nu¬ 
merically by considering small regions over which g e s(t) 
is approximately constant. For large mixing times it can 
be assumed that ho = 0. 

The cocoon surface can only become Rayleigh-Taylor 
unstable some time after the corresponding ensemble el¬ 
ement has entered the subsonic expansion phase. This 
ensemble element will thereafter entrain material from 
the external medium in the mixing region. The thickness 
of the mixing region will grow along the radial axi^l, as 
described by Equation IIII The mixing layer is assumed 
to move with the working surface, expanding equally in¬ 
ward and outward from this equilibrium radius. The 
visible radius of the cocoon R v at a given angle 9 at 
each time step is thus R v (9) = R(9 ) — \h{t, 9). The 
visible, unmixed volume of the cocoon element is thus 
dVy(9) = ^ttR v 3 (9) sin 9d9 (c.f. Equation [3]) . The frac¬ 
tion of this total cocoon volume that is unmixed and thus 
visible through radio emission is 

t - _ / dV v (9) _ jR v 3 (9)sm9d9 

f dV(9) f R 3 (9) sin 9d6 { ’ 

3 The initial exponential growth phase is insignificant in com- 
parison to the self-similar growth. Using linear stability theory of 
I Rayleigh! d 18831) we found that a delay of 10 e-folding times yielded 
a less than 10% variation in the unmixed volume and hence lumi¬ 
nosity. 

4 The same result is obtained if we instead assume that the 
mixing layer grows normal to the cocoon surface with the effective 
gravity reduced by the angle 6 between the radial axis and the 
normal (i.e. g e ff± = g e f[Cos6). For large times t^> to, the mixing 
layer thickness is proportional to g e ff , and thus the thickness along 
the radial axis is hj_(t)/ cos 0 = K,g e ffj_/ cos# = K,g e ff = h{t). 
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Figure 4. Visible cross-section for a radio source with Qtot = 
1q35.5 w total j e t power and initial axis ratio of A = 2, expanding 
into a hos t galax y with typical external pressure profiles described 
in Section 14.2.31 The contours on this figure correspond to source 
ages of 3, 6, 12, 25, 50 and lOOMyrs. The cocoon becomes sub¬ 
sonic and begins to entrain thermal material due to the Rayleigh- 
Taylor instability, first along the transverse axis (at ~ 0.2Myrs), 
and later along the jet axis (at ~ 15Myrs). This Rayleigh-Taylor 
mixing changes the shape of the visible cocoon, eventually causing 
its length to decrease with time, after ~ 40 Myrs. 

The radio source luminosity calculated in Section [2] 
is reduced by this factor to yield the luminosity from 
the non-mixed, synchrotron emitting part of the cocoon, 
A[/ V = RL„. We note that this equation assumes uniform 
luminosity through the unmixed regions of the cocoon. 
tThis assumption may not strictly be valid in weak FR-I 
type sources where the luminosity is often dominated by 
the core. 

Simulated cross-sections of the visible, unmixed co¬ 
coon of a radio source are shown in Figure [4] Rayleigh- 
Taylor instabilities affect the slow-expanding narrow part 
of the radio cocoon first, giving the cocoon a charac¬ 
teristic hourglass shape observed in FR-II sources such 
as Cygnus A. The Rayleigh-Taylor mixing may also 
explain the observed pinching and eventual separation 
of radio bubbles from their active nucleus (e.g. M87; 
iForman et al.l[2007h . In Section [SJ we explore further the 
relationship between cocoon entrainment and observed 
morphology. 

4. ENVIRONMENT AND PARAMETER ESTIMATION 

In this section, we describe the methodology by which 
we use our transonic dynamical model to constrain the 
physical properties of observed AGN. The intrinsic AGN 
parameters, in particular their jet powers and source 
ages, are determined for each radio source based on their 
observed size and luminosity. We quantify radio source 
environments using semi-analytic galaxy formation and 
evolution models. 

4.1. Local radio-loud AGN sample 

IShabala et all (120081) constructed a complete, volume- 
limited sample of radio AGN at 0.03 ^ z ^ 0.1 by cross¬ 
matching two VLA radio surveys at 1.4 GHz, the NR A O 
VLA Sky Survey (NVSS; iCondon et al.|[T998h and Faint 
Images of the Radi o Sky at Twenty-Centimeters (FIRST; 
iBecker et al.ll!995h , and the optical S PSS Data Re l ease 2 . 
We further pair this sample with the lYang et~ahl (|2007T l 
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galaxy clustering catalogue to obtain an estimate of the 
galaxy stellar mass and the mass of the halo associated 
with the cluster. The masses of the central black holes 
are estimated using the iGiiltekin et al.l (2009) determi¬ 
nation of the Mbh ct* relation for elliptical galaxies, 
where er* is the stellar velocity dispersion. This yields 
a sample of 615 radio-loud AGN with stellar and halo 
masses. _ 

iBest fc Heckmanl (|2012l ) analysed a similar (but larger) 
sample of AGN, and separated these into high- and 
low-excitation radio AGN. Cross-matching our catalogue 
with theirs, we obtain nine HERG and 336 LERG classi¬ 
fications, with the remaining 270 galaxies unclassified. 
Both populations are mostly of elliptical galaxy mor¬ 
phology (HERGs 88% and LERGs 93%), and we hence 
adopt environments of simulated elliptical galaxies here 
(see Section [4.2.111 . The emission line fluxes of the ra¬ 
dio source host galaxi es are taken from the value-added 
MPA-JHU catalogues dTremonti et al.ll2004l) . In particu¬ 
lar, we use the [O III] 5007 line to calculate the accretion 
rate of our AGN. 

4 . 2 . Radio AGN environments 
4 . 2 . 1 . Hot gas density 

Environments into which radio sources expand play an 
important role in determining the observed AGN prop¬ 
erties. X-ray observations are not available for the ma¬ 
jority of our sample, and in any case would be biased 
to rich environments. We therefore turn to the semi- 
analytic galaxy evolution (SAGE) model of lCroton et all 
(l2006f) to estimate realistic values for the total hot gas 
mass based on the mass of the cluster halo or galaxy 
sub-halo. The semi-analytic prescription for the growth 
of galaxies and their central supermassive black holes is 
implemented on t op of the output of t he Millennium N- 
body simulation (|Springel et al.l 120051) . and is available 
from the Theoretical Astronomical Observator In this 
work we use a subset of the mock galaxy catalogue span¬ 
ning a simulation volume 30/i -1 Mpc on the side. The 
mock galaxy catalogue consists of both elliptical and spi¬ 
ral galaxy populations. Radio-loud AGN in our sample 
are preferentially hos ted by massive elliptical galaxies 
CSection Id. 11 see also IBest et al.l f2005h . and so we only 
consider these galaxies in quantifying radio source envi¬ 
ronments. 

We compare the hot gas masses obtained from 
SAGE to observations of the gas mass in clusters 
dMcGaugh et al.1 [20 1 01 iGonzalez et al.ll2013l! . The ob¬ 
served gas masses in these > 10 135 M 0 halo mass clus¬ 
ters are typically a factor of two lower than predicted 
by SAGE. We therefore lower the SAGE gas masses by 
this factor when using our model to estimate the physi¬ 
cal parameters of AGN. In Section POl we investigate the 
sensitivity of our results to this choice of hot gas masses. 

An interesting question is whether for AGN that are 
not hosted by brightest cluster galaxies (BCGs), the rel¬ 
evant hot gas mass is that of the host galaxy sub-halo, 
or of the whole cluster. We fit a typical density profile 
(see Section 14.2.31) to the simulated gas mass and virial 
radius of such haloes. The gas density contribution at 
a typical radio source size of 10 kpc from the active nu- 

5 https://tao.asvo.org.au/tao/ 
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Figure 5. A comparison of main and sub-halo hot gas densities 
for the galaxies in our sample. The distance of these galaxies from 
the centre of their group or cluster (scaled by the virial radius) is in¬ 
dicated by the point colour. The straight line of points corresponds 
to brightest cluster galaxies whose sub-haloes are associated with 
the majority of the gas in their cluster. 

cleus is modelled for both the cluster halo and sub-halo 
using simulated galaxies in SAGE with comparable stel¬ 
lar masses to the bulk of our sample (i.e. > 10 la5 M 0 ). 
These two contributions to the gas density are plotted in 
Figure [3 and grouped by the radial location of the host 
galaxy within the cluster halo, r/r 2 oo- Approximately 
80% of these are BCGs (black points in the figure), but of 
the remainder, less than 2% have environments at 10 kpc 
with a greater gas density contribution from the cluster 
halo profile than the host sub-halo. Subsequently we 
therefore model the external environment using the sim¬ 
ulated hot gas density of the host galaxy sub-halo. 

4 . 2 . 2 . Gas temperature 

The temperature of the hot gas in clusters is known 
to scale with the hal o mass as T oc M 2 oo°' 67± °'° 6 (e.g. 
lYikhlinin et al.ll2006D . in line with the virial tempera¬ 
ture of the halo. However, the gas temperatures in this 
relationship exclude the central cooling cores and do not 
includ e any recent heating from AGN. lO’Sullivan et al.l 
(120011) compiled a catalogue of 401 early-type galaxies, 
including 13 AGN hosts, with X-ray luminosity obser¬ 
vations as a function of the R-band optical luminos¬ 
ity. We c onvert these var iables to temperature and 
halo mass (jBell et all 12003 ) to obtain the dependence 
T oc M 2 oo°' 27±0 ' 24 - If we include only the ten AGN with 
haloes in the mass range of our sample we find no signif¬ 
icant relationship (T oc M 2 oo 0 ' 1S±0 ' 38 )- The temperature 
of the AGN environment may therefore have only a weak 
dependence on the mass of its cluster halo. At the me¬ 
dian halo mass of our sample, the best fit temperature 
based on the X-ray luminosity of the AGN is 10 70± °- 3 K. 
In this paper we present results that assume this constant 
mean gas temperature. We note that our findings are rel¬ 
atively insensitive to the choice of temperature profile, 
since very few of our AGN hosted by low-mass galaxies 
are expected to be in the subsonic expansion phase where 
the effects of ambient pressure may be important. 

4 . 2 . 3 . External pressure profile 

The radial dependence of the host gas density and 
temperature, which together form the overall pressure 






































Table 2 

Ambient pressure profile variables. 


i 


Fixed parameter 

Value 


initial axis ratio 

A 

2 


slope (steepening) 

£ 

3 



i' 

3 


steepening radius 

r s /r200 

0.7 


Sampled parameter 

Mean 

(7 

slope (core) 

a 

1.7 

0.3 

core radius 

ao/r200 

0.08 

0.03 

slope 

P 

0.6 

0.1 

temperature 

loglO r x 

7.0 

0.3 

gas mass 

Afgas 

mock galaxy data 

virial radius 

^200 




profile, is modelled using the cluster observations of 
IVikhlinin et al.l (1200611 . Their sample comprises low red- 
shift (z < 0.3) clusters with halo masses around 10 14 M 0 . 
The density profiles of these clusters are observed to have 
very similar shapes, when scaled to the core density and 
virial radius of the cluster. The temperature of the clus¬ 
ters is similarly constrained, but show far less variabil¬ 
ity over the radius of the cluster than the density (i.e. 
|£| <C |/31, using parameters in Table [[]). The mean gas 
temperature (Section 14.2.21) is therefore adopted for all 
radii, with £ taken as zero. 

The gas density profile assumed in this work is 
a simplified fo r m of the double-/3 profile derived by 
IVikhlinin et al.l (120061 their Equation 3) for a sample of 
13 local clusters. Only six of their thirteen clusters need 
a second, core region /3 profile in their fit, and as such 
we have chosen to neglect this term here to reduce the 
degrees of freedom. Our simplified hot gas density profile 
is thus 


(r/a 0 ) a _ 1 _ 

(1 + r 2 /a 0 2 ) 3/3 '““/ 2 (1 + r"<'/rj'yh' 

(13) 

where po is the mass density in the core, a o is the core 
radius and r s > 0 . 3 r 2 oo is the radius at which the profile 
steepens. The exponents a, /?', 7 ' and e set the slope 
of the profile, whilst ao and r s are the core radius and 
the radius at which the double-/? profile steepens. The 
means and standard deviations we assume for these vari¬ 
ables in this work are summarised in Table [2] In our 
radio source modelling, this continuous density profile 
is approximated by a series of ~ 100 power laws (see 
Equation [1]), conserving the total gas mass M gas of the 
sub-halo within r 2 oo as simulated using SAGE for a given 
host stellar mass. In Section 14.3.21 we show that the de¬ 
rived physical parameters for observed radio sources are 
relatively insensitive to realistic changes in these external 
pressure profiles. 

4.3. Jet powers and ages 

The external pressure profile p x (r) of each observed 
galax y is model l ed by a set of 128 profiles based on the 
IVikhlinin et al.l (120061) density and temperature profiles 
with parameters given in Table [2] For each profile, we 
sample the gas mass and virial radius distributions from 





Jet power [log W] Source age [log yrs] 


Figure 6. Marginalised, normalised one-dimensional probability 
distributions for the jet power (left) and source age (right) for 
typical AGN in our sample. 


the simulated SAGE galaxies, at the stellar mass of in¬ 
terest, and normal distributions (shown in Table O on 
parameters describing the shape of the density and tem¬ 
perature profiles. We simulate radio source evolution for 
a range of jet powers Q tot and ages f age in each pressure 
profile. The goodness of fit of the simulated luminosity- 
size profile (or L- D track) for each such environment 
to the observed galaxy size D a b s and luminosity T 0 bs (in 
log-space) is assessed using the y 2 statistic, 


x 2 = 


-^pred -^obs A f -^pred -^obs \ 




■ 


0L 


(14) 


where _D pre d and L pre d are the corresponding parame¬ 
ters predicted by our model. Here the uncertainty due 
to the model dominates and is taken as < 7 n = ay = 
0.3 dex. The likelihood of each simulation matching the 
observations is calculated as P pre d = e _x / 2 , yielding 
a joint probability distribution P(Qtot,t age: p x ). The 
marginalised probability for the n-th value of the sim¬ 
ulated jet power Qtotn matching the observations is then 


P(Qtotn) — EE P(Qtot n? ^age i? Pxj), (15) 

i =1 3=1 

with an equivalent expression for the source age, 

-^(^agen) — EE P(,Qtot i? ^age n? Pxj). (16) 

i=1 3=1 

These probability distributions are approximately log¬ 
normal (Figure [ 6 j). The best estimates for the jet power 
and source age are thus taken as the peaks of the relevant 
marginal distributions, with the uncertainties taken as 
the FWHM. We repeat this procedure for each of the 
615 AGN in our sample. 


4.3.1. Robustness to jet power and age priors 

The joint probability distribution is initially calculated 
assuming a flat prior (in log-space) for the jet power and 
source age, and a log-normal distribution for the host 
gas density. The estimated distributions of these param¬ 
eters are coupled (e.g. both old, high jet power, and 
young, low jet power sources will appear faint and are 
thus more likely to be undetectable in our radio surveys). 
We therefore relaxed the assumption of log-uniform pri¬ 
ors, instead adopting the estimated distributions (from 
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the first pass) of jet powers and ages in the marginalisa¬ 
tion process. The results are robust to these changes in 
the priors on jet power and source age. 

4.3.2. Robustness to environment 

We further examine the reliability of our jet power and 
source age estimates by assessing their sensitivity to the 
host gas density prior. The jet powers and source ages 
are recalculated for gas density priors with a factor of two 
shift (i.e. 0.3 dex) in the mean density of the distribu¬ 
tion at each stellar mass. The shifts in the estimated jet 
power and source age distributions are shown in Figures 
I10l through ll3l The higher pressure prior yields lower jet 
powers, in particular for our lower mass hosts, and older 
sources in all galaxies. This is expected sin ce h i gher ex¬ 
terna l pressures yield higher luminosities (jKaiser et alJ 
119971 Equation 16) requiring a lower jet power to match 
observations. The source age is necessarily greater due 
to the increased resistance to expansion. Conversely, the 
lower pressure prior corresponds to higher jet powers and 
younger sources. These factor of two shifts in the mean 
of the external pressure distribution result in a less than 
0.15 dex variation in the median jet power and source 
age. These variations are both much less than the in¬ 
trinsic width in the derived jet power and source age dis¬ 
tributions (0.6 and 0.9dex respectively). We conclude 
that the jet power and source age parameter estimates 
are therefore robust to the moderate uncertainties in the 
distribution of the pressure profiles. We estimate these 
parameters are accurate to better than a factor of 1.5 
(in a statistical sense). We note that our assumption of 
a cluster-like pressure profile may not be applicable for 
low-mass AGN hosts. However, the assumed gas masses 
(from SAGE) are consistent with observed galaxy prop¬ 
erties in a statistical sense, validating our approach. 

4.4. Axis ratio 

The initial supersonic phase value of the axis ratio, A 
(this is the ratio of lobe length and radius), is a free 
parameter in our model. We adopt a value of A = 2 
here for all sources, noting that A increases once the 
cocoon enters the transonic phase. Our model predicts 
axis ratios in the range 2-8, bro adly consistent wit h the 
observed radio AGN population (jMullin ct al.ll2008l) . We 
defer a detailed investigation of cocoon shape to a future 
paper. 

5. RADIO SOURCE MORPHOLOGY 

Thc lFanaroff fe Rilevl morphological dichotomy can be 
at least partly explained by the influence of external en¬ 
vironmental factors on the jet structure (e.g.lLaing et all 
11994 iBicknelll 119951 : iKaiser fc Alexander! 1 199 7th In nar- 
ticular, shredding of the subsonic cocoon by Rayleigh- 
Taylor mixing of the external medium makes the jet 
more liable to Kelvin-Helmholtz instabilities, leading to 
an FR-I like morphology. Here, we test whether the co¬ 
coon expansion speed predicted by our dynamical model 
is related to radio source morphology. In particular, we 
hypothesise that sources expanding supersonically will 
appear as FR-IIs, and those expanding subsonically as 
FR-Is. 

In our model, morphology is determined by calculating 
the fraction of their lives that radio sources of a given jet 



Figure 7. FR radio dichotomy as a function of i?-band opti¬ 
cal magnitude and 1.4 GHz lumino sity. Obser vations are from 
I Owen fc Ledlowl d 19941 , points) and I Best! (120091 . contours). The 
contours correspond to locations with 16, 50 and 84% fewer sources 
than the peak phase space density. These observations are scaled 
assuming our modelled results used f = 1. Higher values in this 
uncertainty factor will cause these points to move to higher lumi¬ 
nosities relative to our predictions. The straight lines correspond 
to model predictions where half the sources are still expanding 
supersonically, and the other half subsonically. Calculations for 
expansion speed along the jet axis (dashed line) and transverse to 
the jet (solid) are shown, assuming a constant 10 7 K temperature. 
The lcr uncertainties due to variations in the temperature of the 
external medium are shown by the grey shading. 

power and stellar mass spend in the subsonic and super¬ 
sonic expansion phases. The source active lifetimes are 
constrained by fitting a multivariate regression to the jet 
power and stellar mass estimates we obtained for our 
sample. Because the cocoon is non-spherical, parts of 
the surface become subsonic at different times. Hence 
the choice of the reference point for measuring the co¬ 
coon expansion speed is important. Here we consider 
two locations, along the semi-minor (or transverse) and 
semi-major (or jet) axes of the cocoon. These are parts 
of the cocoon that become subsonic first (semi-minor) 
and last (semi-major axis). 

We compare this simulated FR-I/II dich otomy with 
two se ts of observations in Figure [7] IQwen fe Ledlowl 
(|1994f l reported the observed correlation between R-band 
optical magnitude, 1.4 GHz radio luminosity and FR 
morphology for a s ampl e of rad io sources primarily found 
in Abell clusters. iBestl (l2009f) pointed out that, in this 
sample, the FR-Is are located predominantly at redshifts 
less than 0.1, with the bulk of the FR-II population be¬ 
ing found at z > 0.25. By constructing a volume-limited 
sample. IBestl (I2009T) showed that the separation between 
FR-Is and FR-IIs in optical magnitude-radio luminos¬ 
ity phase space, albeit real, is much le s s clea rcut than 
originally suggested bv lOwen fc Ledlowl (119941 . 

In Figure El the dashed and solid black lines represent 
fits from our model to the locations in jet power-stellar 
mass space where an equal fraction of the radio source 
active lifetime is spent in the subsonic (i.e. FR-I) and su¬ 
personic (i.e. FR-II) states. These fits are for a model us¬ 
ing a constant 10 7 K temperature for environments host¬ 
ing our AGN, whilst the grey shading indicates the lcr 
uncertainties corresponding to a factor of two variation 
in the temperature of the external medium. We convert 
stellar mass to R-band optical magnitude using the equa- 
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Figure 8. FR radio morphology as a function of source size, 
1.4 GHz lobe luminosity and host galaxy stellar mass. The shad¬ 
ing maps regions in this phase space where our model predicts 
transonic phase evolution. Sources above these dividing lines for 
the relevant mass bin are predicted to be FR-IIs, while those below 
should be FR-Is. Objects lying within the shaded transonic regions 
at a given stellar mass can have either morphology. The morphol¬ 
ogy of sources is not expected to change further upon entering the 
region beyond the grey dashed line. That is, a source retains its 
FR-II morphology if its L-D track does not cross the morpholog¬ 
ical transition line before entering this region. Note the overlap 
between shaded regions for different h ost galaxy masses. Points 
show data from the CoNFIG sample of lGendre et al.l (120101) , with 
FR-IIs marked as circles and FR-Is as diamonds. Colours corre¬ 
spond to the stellar mass of the AGN host galaxy. These observed 
points have been scaled using / = 1. 

tions of lBell et al.l (120031) . and jet power to radio luminos¬ 
ity using our data. Based on our hypothesis, we expect 
sources lying above the semi-minor axis sonic transition 
line to be predominantly FR-IIs, with FR-Is lying below 
the jet axis line. The empirical separation between the 
observed FR-I and -II populations is consistent with the 
lines derived for the two subsonic points, and corresponds 
to a surface location somewhere between these extreme 
cases. Interestingl y, the region where these two popula¬ 
tions overlap in the lBestl ( 2009 ) sample is coincident with 
the parameter space in which our predicted model mor¬ 
phologies are strongly environment-dependent. These 
conclusions are robust to realistic uncertainties in the jet 
power-luminosity relationship, e.g. / < 5. We conclude 
that a comparison of the expansion speed of the cocoon 
to the local sound speed can thus successfully explain the 
observed FR-I/II dichotomy. 

Based on our findings, we tentatively identify radio 
sources as FR-Is or FR-IIs based on the Mach number 
modelled at the time of observation. Sources which are 
entirely supersonic are labelled FR-IIs, and those which 
are entirely subsonic are FR-Is. Sources which fall be¬ 
tween these two categories could have either morphology, 
with the probability of being an FR-I increasing with the 
fraction of the cocoon surface which is subsonic. The FR- 
II population is predicted to consist of both smaller and 
younger sources than the FR-I population, irrespective 
of which tem perature model is used, consistent with ob¬ 
servations bv lBest et al.l (|2007D . 

The locations in luminosity-size space in which FR- 
I and -IIs are thus expected to be found are plotted in 
Figure [8] for four host galaxy stellar masses. For each 
mass bin, the shading marks the regions in luminosity- 
size at which the expansion of radio sources is in the 


transonic phase (i.e. supersonic at semi-major but sub¬ 
sonic at semi-minor axis). The shaded regions include the 
uncertainty due to simulated variation in the host envi¬ 
ronment within each mass bin. At a given mass and size, 
radio sources in the supersonic expansion phase (which 
we identify with FR-IIs) have higher luminosities than 
these dividing lines, whilst sources expanding subsoni- 
cally (our FR-Is) are less luminous. Sources lying within 
the shaded regions defined by the semi-major and semi¬ 
minor axis sonic transitions can have either FR morphol¬ 
ogy in our model. Radio AGN inhabiting poor envi¬ 
ronments are more likely to be observed as FR-IIs than 
sources of the same size and luminosity in more mas¬ 
sive hosts. However, the evolution of the larger sources 
is dictated mostly by the higher density environment at 
smaller radii. Sources which have reached sizes beyond 
the grey dotted line in Figure [8] without entering the 
transonic phase are thus expected to retain their FR-II 
morphology. The morphology of the largest radio AGN 
is therefore somewhat sensitive to the initial conditions, 
in particular their host environment. We compare our ra¬ 
dio morphology predictions with obs ervations of the loca l 
extended radio source populations (jGcndre et al.112010 ). 
and find the two to be in good agreement. That is, all 
their FR-IIs are in the supersonic or transonic expan¬ 
sion phase, whilst FR-Is are in the transonic or subsonic 
phase. This result holds true for realistic values of the 
scaling factor / in the jet power-luminosity relationship. 

Finally, we note that our data can be used to derive 
an empirical relationship for the jet power based on the 
source morphology, size and luminosity observables. A 
number of loss processes associated with the synchrotron 
emitting electron populations render luminosity (by it¬ 
self) a useful but flawed pr edictor of j et power espe cially 
for evolved radio galaxies dShabala k, Godfrevll2013D . We 
defer a detailed investigation of the jet power-radio lu¬ 
minosity relation to a future paper, but now examine its 
dependence on the Fanaroff-Riley morphology. The jet 
power (calculated in Section H~3l) is plotted in Figure [9] 
as a function of size, luminosity and our tentative source 
morpholog y classification. A l so plo tted is the cavity rela¬ 
tionship of lHcckman fc Best ( 2014 ), scaled assuming our 
modelled results had a jet power-luminosity uncertainty 
factor of / = 5. That is, our model is in good agree¬ 
ment with other reported measurements of the jet power 
when assu ming f = 5 . This value agrees well with the 
findings of iDalv et al.l ( 2012 ) using a different method 
of estimating jet pow er, and is well within th e plausible 
range / ~ 1 - 20 (e.g. iHeckman fe Bestll2014l) . 

The jet power in Figure [9] shows a strong positive cor¬ 
relation with the 1.4 GHz radio luminosity by itself, i.e. 
brighter sources tend to have higher jet powers. More 
of the variability in the jet power can be explained by 
considering the source size and morphology. The source 
size is plotted in Figure |H] for two size bins whilst the 
morphology is indicated by the point colour, with FR- 
Is in red and FR-IIs in blue. Larger sources tend to 
have higher jet powers than smaller sources of the same 
luminos ity and mor p hology, consis tent with the observa¬ 
tions of lShabala fc Godfrey! (|2013t ). Further, radio AGN 
of FR-II morphology have jet powers approximately a 
factor of two higher than FR-Is of the same size and lu¬ 
minosity (dashed lines in Figure[9]). This is likely because 
at a given jet power, AGN expanding in poor environ- 
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Figure 9. Jet power as a function of the source size, luminos¬ 
ity and predicted Fanaroff-Riley morphology. Radio AGN with 
sizes D < 10 kpc a marked by circles, whilst the largest sources 
(D > 10 kpc) are shown with crosses. FR-Is are plotted with red, 
unfilled points and FR-IIs using blue, solid points. Here, FR-IIs 
are identified with sources that are still expanding supersonically, 
while FR-Is are those sources for which at least part of the cocoon 
is expanding subsonically. Best fits to the smallest sources for both 
morpholo gies are p l otted wi th das hed lines. The cavities relation 
derived bv IHeckman fc Best! 11201411 is shown for comparison (grey), 
scaled using / = 5. 

ments are more likely to retain FR-II morphology; these 
sources will also be less luminous than their disrupted 
FR-I counterparts located in denser environments. Dif¬ 
ferent empirical relationships will therefore apply for the 
two FR morphological classes, and thus radio AGN mor¬ 
phology must be considered in addition to size and lumi¬ 
nosity in estimating jet power. 

6 . AGN ENERGETICS AND FEEDBACK 
6 .1. Jet powers and lifetimes 

Having estimated AGN kinetic powers and lifetimes 
(Section rOll . we can now quantify the energetics of AGN 
feedback on galaxies. When drawing conclusions about 
the overall AGN population from our volume-limited 
sample we must be careful of selection effects though. 
Low jet power sources (Q to t < 10 34 W/ 3 / 2 ) are theoret¬ 
ically undetectable by the FIRST/NVSS surveys at any 
time in their evolution, whilst higher jet power sources 
will not be detectable at later times when loss mecha¬ 
nisms (in particular the inverse Compton scattering of 
CMB photons) and Rayleigh-Taylor mixing cause the 
luminosity to fall rapidly. We quantify the fraction of 
undetected sources by modelling AGN with a range of 
jet powers, ages and environments. In our modelling, we 
assume that the AGN environments are well described 
by the scaled SAGE hot gas masses (Section 14.21) , and 
we adopt the observed mass dependence of the radio- 
loud AGN fraction to quantify AGN intermittency. The 
figures in the following sections indicate this drop-out re¬ 
gion with blue lines corresponding to detectable fractions 
of 84, 50 and 16% of radio sources. 

6.1.1. Jet powers 

A combination of theoretical and empirical considera¬ 
tions suggest that the AGN jet power should to be largely 
independent of the host galaxy stellar mass. Jet genera¬ 
tion models predict that the jet power increases approxi¬ 
mately with the black hole accretion rate Mbh oc toMbh 


Figure 10. The correlation between jet power and observed stel¬ 
lar mass. Best fit is shown by a solid line, and uncertainties due 
to a factor of two variation i n hot gas mass by th e shaded region. 
The relationship obtained bv IShabala et ahl (120081) using their sim¬ 
plified model for these same radio galaxies is shown in purple. The 
blue lines show the region of jet power-mass space where 16, 50 
and 84% of sources (from top to bottom) are expected to have 
luminosities below the survey detection limits. 

(|Meiedl200lD . The general consensus is that in massive 
ellipticals at low redshift, the active nucleus feeds en¬ 
ergy back into its host galaxy at a ra te balancing the 
loss of energy through radia t ive coo ling (Best et alJl2007t 
iMcNamara fe Nulsenl 120071 : iFabianl 120121) . Naively, we 
might expect this instantaneous accretion rate to balance 
the gas cooling rate, which scales with black hole mass 
(and thus stellar mass) as M coo \ oc Mbh 1 ' 5 dBest et al.l 
120051) . The jet power would thus be expected to be pro¬ 
portional to the black hole and stellar masses as Q to t oc 
Mbh 1 ' 5 oc M* 1 ' 5 . However, this jet power-stellar mass 
relationship needs to be modified for the AGN feedback 
cycle. Since there may be multiple AGN events within 
a single cooling time, we consider the time-averaged fu¬ 
elling rate for the active nucleus, /rlMbh, where /rl is 
the radio-loud fraction. If the AGN are fuelled by the 
gravitational instability, /rl is e xpected to scale with 
black hole mass as /rl oc Mbh 1 ' 5 dPope et al.l 120121 ). In 
heating-cooling equilibrium (i.e. /rlMbh oc M coo \) the 
jet power should therefore be related to the black hole 
and stellar masses as Q to t ~ M coo i//rl oc Mbh° oc M*°. 

The jet powers calculated for our sample are plot¬ 
ted as a function of stellar mass in Figure [TUI Unsur¬ 
prisingly, a wide range of jet powers are observed at 
a given stellar mass. The relationship between these 
variables is analysed with a standard least squares lin¬ 
ear regression (black line) yielding Q to t oc M*°' 1±al . 
We also fit an LI (least absolute errors) linear regres¬ 
sion to provide a more robust estimate of the slope of 
Qtot oc M+° ,0=l=0 1 . By contra s t, usi ng their simpler FR- 
II only model IShabala et alJ (|2008D found a steeper re¬ 
lationship of Qtot oc M*°' 6±q ' 2 , but with a comparable 
mean jet power (purple line). Further, in this figure we 
show the region of jet power-mass space where sources 
are expected to have luminosities below the survey detec¬ 
tion limits (blue lines and shading). Although detection 
limits may contribute to the lack of mass dependence 
for low jet powers (Qtot < 10 34 5 W / 3 / 2 ), the top end of 
the distribution corresponds to the real AGN population. 
This relation shows no dependence on the host galaxy 
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mass, in agreement with the theoretical expectation as¬ 
suming long-term heating-cooling balance. Our results 
are therefore consistent with the maintenance nature of 
AGN feedback in the local Universe. 

6 .1.2. Active and quiescent timescales 

In Section [4] we derived the age of the radio source at 
the time of observation, not the total time over which 
the source is active. As sources age, they become less 
luminous due to a combination of synchrotron, inverse 
Compton and adiabatic losses. Because source size also 
increases with age, the surface brightness of the oldest 
sources drops rapidly, taking many of these close to (or 
below) the survey sensitivity threshold. Through the 
use of our marginalisation we may select very long-lived 
sources early in their active lifetime whilst they are de¬ 
tectable. If these sources were active up to their max¬ 
imum detectable age we would expect the source ages 
to follow a uniform, linear distribution. However, more 
than 90% of our sources are observed to be younger than 
half their maximum detectable age, with more than 70% 
at a tenth of that age. There is therefore a real over¬ 
abundance of young sources. We conclude that these 
young sources, and all other sources, must be observed 
somewhat close to their active age. The source active 
lifetime is thus taken to be double the observed age, i.e. 
ton = 2f age , since the average source will be observed 
halfway through its evolution. 

The relationship between the derived active lifetime 
and stellar mass for our radio AGN sample is shown 
in Figure [TT| The mass dependence of the active life¬ 
time is again calculated using a linear regression yielding 
t on = (4.8 ± O.5Myr/- 1 / 2 )(M*/lO 11 M 0 )°- 7±01 . The 
derived normalisation for our sample is higher than the 
lMyr/ -1 / 2 lifetime estimated bv IShabala et all (|2008l ) 
for a 10 11 M 0 host galaxy, though the slope is compara¬ 
ble given their large uncertainties. Their younger ages 
result from the use of a simpler strong-shock limit su¬ 
personic expansion model for this mostly FR-I popula¬ 
tion. Finally, the AGN active lifetimes should scale as 
ton = /RLthost/-N, where thost is the total age of the host 
galaxy group or cluster an d N is the numb er of AGN 
outbursts during this time dPope et al.ll2012D . The trig¬ 
gering frequency is thus related to the stellar mass as 
Wtrig = N /thost = /rl/ l on oc In general, the 

radio AGN in the most massive galaxies are therefore 
triggered more frequently, and are also active for longer 
in a given epoch. 

The quiescent timescales (i.e. the time the AGN is 
inactive) are also relevant to understanding the AGN fu¬ 
elling and triggering mechanisms. The active lifetime 
can be used to determine the quiescent timescale as 
iquiet = ion(l - /rl)//rl- Here we use the s caled rela¬ 
tionship for the duty cycle of lBest et al.l (1200511 . The qui¬ 
escent timescales (plotted in Figure fl^l) exhibit a strong 
negative correlation with stellar mass, driven largely by 
the strong dependence of /rl on M*. Many of the radio 
sources hosted by galaxies with stellar masses less than 
M* « 10 11 M 0 are expected to be inactive for longer 
than the Hubble time. A large fraction of these low 
mass galaxies may therefore have never been in the ac¬ 
tive state. We note that here we used a source specific 
calculation of the active age, but a population average for 



Figure 11. The active lifetime as a function o f ste llar mass. The 
blue lines are for the same fractions as in Figure ITU1 but with limits 
for both young and old sources. 



Figure 12. The distribution of quiesce nt t imescales plotted as 
a function of stellar mass, as for Figure 1 11 1 The black dashed 
horizontal line is the Hubble time. 

the quiescent phase, which may not be applicable for in¬ 
dividual sources. We find the quiescent time scales with 
mass as f qu i e t oc M* _0 ' 9±01 . In the absence of a heat¬ 
ing source, we wo uld expect f qu iet oc M coo i _1 oc M* -1 ' 5 
(|Best et al.l 1200511 . The fact that the observed relation¬ 
ship is not as steep may suggest that feedback from pow¬ 
erful radio sources in massive galaxies plays an important 
role, flattenin g out the relat ionship at high masses as first 
proposed bv IShabala et al.l (120081) . 

6.1.3. AGN energetics and feedback 

Examining the energetics of AGN reveals which sources 
provide the most feedback in the low-redshift universe. 
The energy input by an AGN into its surroundings over 
an active phase is equal to the product of its jet power 
and the duration of the AGN event, E to t = QtoUon- The 
input energies thus calculated for our sample are plotted 
in Figure 1151 as a function of stellar mass. The relation¬ 
ship between these variables is again analysed using a 
linear regression yielding E tot oc M*°' 8±0 ' 2 . The energy 
injected by an AGN into its host in an outburst thus 
increases approximately linearly with the galaxy mass. 
Since the AGN in massive hosts are also triggered more 
frequently, we find the average energy injected by the 
central nucleus over time to scale even more steeply with 
mass as E tot oc w t rig-Etot oc M* 1,6±0 ' 2 . This scaling is 
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Figure 13. The energy input of the AGN into their surroundings 
as a function of stellar mass. Detection limits are not shown since 
the effect of selection bias is similar for all masses, with approxi¬ 
mately 30% of sources undetected everywhere. While at a given 
stellar mass we observe a broad distribution in the total AGN en¬ 
ergy output, there is a robust mass dependence for the average 
AGN. 

similar to the theoretical expectation for a cluster in 
heating-cooling equilibrium of E tot oc M c00 i oc M* 1 ' 5 . 
The implications of these findings are discussed in more 
detail in Section ECU Massive hosts thus provide more 
feedback to their surroundings than their lower mass 
counterparts, as required by galaxy formation models 
(|Croton et al.ll2006t iBower et alJl2006lh 

We investigate which AGN provide the most feedback 
in the low-redshift universe by considering the number 
density of sources at a given stellar and radio luminosity, 
and their time averaged energy input, E tot . The sources 
are binned in stellar mass and radio luminosity to yield 
luminosity functions (including V/V max corrections) at a 
fixed stellar mass. The number densities in each bin are 
then weighted by the energy input of each radio AGN 
on a source-by-source basis. This energy flux (i.e. en¬ 
ergy input per unit volume) is plotted as a function of 
radio luminosity in Figure [TI] for two stellar mass bins 
and our entire sample. We plot the energy flux includ¬ 
ing corrections for the survey drop-out at low jet powers 
(solid lines) in addition to the uncorrected data (points 
and dashed lines). Poisson uncertainties (la level) are 
plotted for the raw data. 

The energy flux of radio AGN in high mass galax¬ 
ies (M* > 10 114 Mq) increases strongly with luminos¬ 
ity until Li .4 = 10 25 WHz _1 where it peaks. These 
high luminosity sources contribute a factor of ten times 
more energy to the total galaxy population than those 
at 10 23 WHz _1 . By contrast, the energy flux of radio 
sources in lower mass hosts has only a very weak de¬ 
pendence on luminosity. These sources contribute a fac¬ 
tor of four less energy to the low-redshift galaxy pop¬ 
ulation than those in the most massive hosts, despite 
their greater number density. These result hold when 
corrections for the survey selection effects are applied. 
The lower energy flux associated with low luminosity 
radio sources and those in less massive hosts implies 
the amount of feedback energy such sources provide on 
the total galaxy population is much less than for the 
most massive galaxies near the ~ 10 25 WHz _1 luminos¬ 
ity peak. Observing radio sources at even lower luminosi¬ 
ties may therefore be unnecessary in further understand¬ 



Figure 14. The radio AGN energy flux (energy input by AGN 
into their hosts per unit volume) as a function of the observed 
1.4 GHz luminosity for two mass bins and our entire sample. We 
plot both the raw data (points and dashed lines) and also include 
curves that correct our data for survey selection effects (solid lines). 

ing the effect of AGN feedback on their host galaxies. Of 
course this does not take into account the very likely pos¬ 
sibility of different coupling efficiencies with the gas for 
supersonic and subsonic cocoons. As discussed in Sec¬ 
tion [5j lower mass galaxies preferentially host supersonic 
phase FR-IIs whilst their more massive counterparts har¬ 
bour subsonic FR-Is. Even for a given radio source mor¬ 
phology, numerical simulations of radio jets suggest the 
efficiency of feedback is strongly dependent on how the 
energy is suppl ied to the ga s, rather than just the total 
amount dBinnev et al.ll2007l) . 

6 .2. AGN heating and cluster cooling 

In this section, we compare the time-averaged rate of 
heating from the AGN with the cooling in its host’s hot 
environment. Our jet powers calculated on a source-by¬ 
source basis provide a much broader distribution at a 
fixed stellar mass than previous wor k using em pirical 
jet power-luminosity relationships fe.g. IBest et al.|[2007l f . 
We estima t e the cooling rate, following the method of 
iBest et all (|2007lh by considering the fraction of X-ray 
luminosity arising from within the cooling radius, r coo i- 
This rate is found as a function of the clus ter veloci ty 
disper sion, cr v , by combinin g the results of iPeres et al.1 
(1199811 and IQrtiz-Gil et al.l (12004T) . L x (r < r coo i) ss 
1.34 x 10 37 (< 7 v / 1000 kms -1 ) 4 - 1 W. Here, the cluster ve¬ 
locity dispersion is derived from the halo mass of eac h 
cluster using the relationship of iMunari et al.1 1 2012 ). 
The average rate of heating from the AGN is calculated 
using our results as the product of the jet power and 
the observed radio-loud fraction, i.e. H = Qtot/iu/Af*). 
Ag ain we assume the black hole-duty cycle relationship of 
IBest et all (I2005H. but tak e /rl = 1 for the most massive 
black holes lhh,bianll2012l and references therein). Based 
on our findings in Figure O we assume the uncertainty 
scaling factor in the jet power luminosity relationship is 
/ ~ 5 when comparing the heating and cooling rates. 

The ratio of the time-averaged heating to the cluster 
cooling rate is plotted in Figure [15] Low mass galaxy 
haloes are consistent with a ratio of approximately unity 
(i.e. heating-cooling balance) for / < 5, however, high 
mass haloes (> 10 14 5 M 0 ) require an additional heating 
mechanism. Such clusters have time-averaged heating 
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Figure 15. The ratio of the AGN heating and gas cooling rates 
in the host cluster as a function of halo mass. The orange circles 
are the ratio of the time-averaged energy input from the jet and 
the cluster cooling rate. The solid black line is a best fit to this 
distribution with the grey shading showing uncertainties arising 
from our assumptions about the hot gas density. Heating due to 
thermal conduction is plotted in red for two scaling factors. 

rates a factor of 200/ _3 / 2 lower than the cluster cooling 
rates, though due to the large scatter in jet powers this 
difference is only significant at the 1 <t level for f ~ 5. 
This d iscrepancy was previously reported by iBest et SI 
(l2007lf . 

Thermal conduction of energy from the reservoir of hot 
gas outside of the cooling radius into the cooler cluster 
centre may provide an additional source of heating in 
the most massive systems (see IBest et al.l[2007[ and ref¬ 
erences therein). The heating rate due to this thermal 
condu ction wa s origi nally calculated for a pure hydrogen 
gas bv ISnitze'i] (jl962H . Although t his rate w ill be reduced 
i n the presence of magnetic fields, iNaravan fc Medvedevl 
(|2001lf argue that the suppression of thermal conduction 
may only be a factor of a few in the presence of turbulent 
magnetic fields. In Figure[l5]we plot the effects of Spitzer 
conductivity for 25 and 100 per cent of the Spitzer value. 
The inclusion of this heating mechanism for realistic val¬ 
ues of Spitzer conductivity can yield heating-cooling bal¬ 
ance (within lcr uncertainties) for all halo masses. That 
is, AGN feedback, in conjunction with thermal conduc¬ 
tion, can maintain long term heating-cooling balance in 
even the most massive gala xy clu sters, consistent with 
the findings of IBest et all (j2007l f using a different jet 
power estimation method. 

7. FUELLING THE MONSTERS 

In this section, we use our results to study black hole 
accretion flow mechanisms. This requires careful consid¬ 
eration of observational selection effects. 

7.1. High- and low-excitation radio galaxies 

First, we investigate whether the different black hole 
fuelling mechanisms of high- and low-excitation radio 
galaxies (HERGs and LERGs) correspond to differ¬ 
ences in the active lifetime and jet power distributions, 
as might be expected if these objects probe different 
black hole accretion s tates (e.g. lHardcastle et all 120071 : 
IBest fe Heckmanll2012lf . There are not enough HERGs in 
our sample to draw reliable conclusions, however, the rel¬ 
ative locations of the HERG and LERG sub-populations 
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Figure 16. The distribution of high- and low-excitation radio 
galaxies, pentagons and circles respectively, in active age-jet power 
space. The black dashed curve is a best fit to the distribution of 
LERGs, used to remove the effect of the active lifetime and se¬ 
lection effects from comparisons with the HERGs. The colouring 
of the points denotes the host galaxy stellar mass, whilst the blue 
lines indicate the times at which lower jet power sources are theo¬ 
retically too faint to be detected (solid-50% detectable, dashed-16 
and 84% detectable). 



Figure IT. Stellar-to-halo mass ratio of HERG (black pentagons) 
and LERG (red circles) populations as a function of halo mass. 
There is no statistically significant difference between the stellar- 
halo mass ratio for the HERGs and the best-fit line (black) for 
the LERG population. The single HERG identified as a sp ir al (i.e . 
SDSS r-band fracDeV parameter < 0.8. 1 Padilla &; Strauss! [20081) 
is circled, and again is consistent with the LERG distribution. 

in active age-jet power space are shown in Figure 1161 
The drop-out region of age-jet power space is marked 
with blue lines corresponding to detectable fractions of 
84, 50 and 16% as discussed in Section lCTTI The observed 
“correlation” is thus largely a result of these survey selec¬ 
tion effects, however the relative locations of HERGs and 
LERGs within this distribution are likely due to intrinsic 
differences. 

The HERG population has a mean stellar mass 
0.24 dex lower than for the LERGs, with this difference 
significant at the 3cr level. This is consistent with HERGs 
being fuelled by minor mergers, and LERGs through gas 
cooling which preferentially takes plac e in the most mas¬ 
sive galaxies (lHa rdcastle et all 120071 IBest fe Heckmanl 
1M1 IShabala et al. 1 1201211 . The HERG and LERG ra¬ 
dio sources, however, show no preference for inhabiting 
different environments for a given halo mass (Figure [T7]). 
Hence, we are not just seeing young HERG AGN because 
they are in poor environments (at a given M*). 
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The mean jet power of the HERG population is higher 
than that of the LERGs by 0.58dex (3er significance), 
as shown in Figure 1161 This result holds even when 
the nine HERGs are compared with a sub-population 
of LERGs with the same distribution of stellar masses. 
These higher HERG jet powers are expected due to the 
higher black hole accretion rates in minor mergers. More¬ 
over, the mean active lifetime of the HERGs is 1.1 dex 
lower than for the low-excitation population, with this 
difference significant at the 5a level. These younger 
HERG active ages likely result from these galaxies be¬ 
ing associated with short lived epochs of high accretion 
arisi ng from the merger of the host with a gas ri ch sys¬ 
tem (lHardcastle et al.ll2007t IShabala et al. 1 1201211 . The 
timescales for the supply of cold gas in this scenario dif¬ 
fer from the usual heating and cooling cycle of the IGM, 
which can fuel the lower accretion rate LERGs. 


7.2. Accretion rates and jet generation efficiencies 

The different fuelling modes of HERGs and LERGs 
are expected to yield different accretion rates. We inves¬ 
tigate the distribution of accretion rates for these pop¬ 
ulations by considering the kinetic power of the active 
nucleus relative to its Eddington luminosity, Lecm- The 
central black hole radiates energy directly from its accre¬ 
tion disk with the bolometric luminosity of the disk L ra d 
in addition to mechanical output of energy through the 
jets Lmech = Q tot- The power emitted by the AGN is 
thus related to the black hole mass accretion rate Mbh 
as (L ra d + T mec h) = cMbhc 2 , where e is the efficiency at 
which the accreted rest-mass energy is converted. The 
Eddington-scaled accretion rate is therefore 


Mbh _ Trad + Tmech 
il^Edd -^Edd 


(17) 


We note that this accretion rate is unaffected by re¬ 
alistic values of the jet power luminosity scaling factor 
/ < 5 since the mechanical jet power is much less than 
the disk bolometric luminosity. The bolometric radia¬ 
tive luminosity of each radio source is estimated using 
the observed lumino sity of the [OIIIl 5007 emission, and 
Trad ~ 3500 Loiii (lHeckman et al.1 200411 . The uncer¬ 
tainty on individual estimates of the bolometric radia¬ 
tive luminosity is ~ 0.4 dex. Although 43% of all radio 
sources in our sample have emission below the formal 
line flux uncertainties at the 3a level, this fraction falls 
to less tha n 10% for objects c l assifie d as either HERGs or 
LERGs bv iBest fe Heckmanl (|2012f) . Following on from 
their work we adopt the 3er upper limit to [O III] lumi¬ 
nosity for such sources. 

In order to compare our accretion rate estimates with 
the theoretical expectations we need to correct for the 
survey selection effects. There are two factors we must 
consider: first, the drop-out of lower luminosity sources 
due to Rayleigh-Taylor mixing and late time loss mech¬ 
anisms; and second, the greater likelihood of detecting 
long-lived sources from the total population compared 
to their younger counterparts. Specifically, the number 
density of sources in the population is probabilistically 
expected to be inversely proportional to its active age, 
rc agn (T>n) oc 1 /ton- The number of radio AGN missing 
due to the drop-out of low-luminosity sources is insignif¬ 
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Figure 18. Top panel: distribution of Eddington-scaled accretion 
rates for LERGs (red) and HERGs (blue). The background shad¬ 
ing shows the theoretically expected values for an ADAF (or slim 
disk; red) and a thin disk accretion flow (blue). The dashed lines 
show the uncorrected observed data, and solid lines are corrected 
for selection effects. Middle panel: mechanical jet power as a frac¬ 
tion of the Eddington luminosity. Bottom panel: efficiency of jet 
production. Fractions of the total AGN population are plotted, 
with HERGs scaled by a factor of ten to aid visibility. 

icant in comparison, as we discuss here. The luminos¬ 
ity of our simulated sources typically peaks at approxi¬ 
mately lOMyrs (see Figure [3j). That is, we are able to 
observe all sources with jet powers above 10 345 W/ 3 / 2 
in this age range, as seen in Figure ITfil By contrast, all 
high jet power sources (> 10 35 5 W/ 3 / 2 ) younger than 
~ 1 Gyr can be detected with our survey. The distribu¬ 
tion of counts of these high jet power sources in active 
lifetime space can thus be scaled to the source count at 
lower jet powers to predict the number of missing ob¬ 
jects. In this manner, we find that at least half of the 
10 34 - 5 < Qtot < 10 35 -° W/ 3 / 2 jet power sources expected 
from the counts of 10 35 5 ^ Q tot < 10 36 ° W/ 3 / 2 sources 
are observed. By contrast, correcting for the likelihood of 
observing objects based on their lifetimes yields a factor 
of ten increase in number density of radio sources for each 
order of magnitude decrease in the active age. Moreover, 
the derived active ages cover a much larger range (4 or¬ 
ders of magnitude) than the jet power (2.5dex). Hence, 
correcting observations for age bias is more important 
than the effects of different jet powers. 

The distribution of accretion rates thus corrected is 
now plotted in the top panel of Figure [lH] for the small 
population of HERGs (blue) and the much larger pop¬ 
ulation of LERGs (red). The sub-population of LERGs 
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has a mean accretion rate of to = 0.0006 with a FWHM 
of 0.8 dex with typical (i.e. > 95% of AGN) values in 
the range 10 -4 < rh < ICG 2 ' 5 . The high-excitation pop¬ 
ulation clearly has higher accretion rates with a mean of 
rh = 0.02, with a difference significant at the 9er level. 
In fact, all but one of these nine objects have accretion 
rates greater than rh = 10 -2 ' 5 compared to less than 2% 
of LERGs. Note that sources with luminosities below the 
flux detection limits of the FIRST and NVSS surveys for 
their entire lifetimes remain absent in these comparisons. 
These accr etion rate distribu t ions a re consistent with the 
analysis of lBest fc Heckmanl (120121) for redshifts z < 0.1 
(shown in their Figure 6), and largely with the theoret¬ 
ical division between the ADAF and thin disk accretion 
flow mechanisms at approximately rh = 0. 01 (red and 
blue background shading in Figure [18] iMeierl [200111 . 

The middle panel of Figure |T8] shows the mechani¬ 
cal jet power scaled by the Eddington rate, Lmech/AEdd- 
The HERGs have Eddington luminosity scaled jet powers 
greater than those of the LERGs by approximately a fac¬ 
tor of twelve (6cr significance). The relatively large frac¬ 
tional difference in accretion rates between the HERG 
and LERG populations compared to their Eddington- 
scaled jet powers indicates that the HERGs are much 
less efficient at producing jets. The jet production effi¬ 
ciency is plotted in the bottom panel of Figure [TB] The 
LERG population has jet production efficiencies typically 
in the range 0.01 < £j e t,LERG// 3 ^ 2 < 0.3 with a mean 
efficiency of ejet.LERG// 3 / 2 = 0.06. By contrast, 82% 
of HERGs have efficiencies less than this value with a 
mean of £j e t,HERG// 3 ^ 2 = 0.02. The high-excitation ra¬ 
dio galaxies are therefore approximately a factor of three 
less efficient at producing jets than the LERGs, and this 
difference is significant at the 4er level. 

7 . 2 . 1 . ADAFs 

The accretion rates and jet production efficiencies of 
the HERG s and L ERG s can be compared to the pre¬ 
dictions of iMeieil (120011 1 for black holes powered by an 
advection dominated or thin disk accretion flow. The 
jet production efficiency of an ADAF powered rotating 
bl ack h ol e is re lated to its jet power Q to t (see Equation 7 
of iMeieil 120011) and the Eddington-scaled accretion rate 
rh by 

e A DAF // 3/2 = QADAF/LEdd = 0.05(0.55+ 1.5j+/). 

m 

(18) 

The bl ack h o le sp in ^ g [0,1] is a difficult qu antity to 
measure. iMartinez-Sansigre fc Rawlings! (120111 ) have es¬ 
timated that in the low-redshift universe the mean black 
hole spin is (j) ~ 0.35, with approximately 24% of black 
holes having spins greater than j = 0.5. Variation in 
this parameter can generate jet production efficiencies 
in the range 0.03 < cadaf// 3 ^ 2 < 0.2, consistent with 
the peak in the drop-out corrected LERG distribution 
(solid red line; FigureQ+l bottom) which we identify with 
ADAFs. The smaller “bump” in the observed distribu¬ 
tion (dashed red line) at an efficiency of approximately 
Cjet.LERG// 3 ^ 2 ^ 0.5 is not consistent with this theoreti¬ 
cal model. These sources have median radio luminosities 
a factor of 45 greater than the rest of the LERG pop¬ 
ulation, whilst sizes, jet powers and active lifetimes are 


25 times larger. The typical sizes of these sources’ radio 
cocoons are of the order of 100 kpc with a median active 
lifetime of 100 Myr/ _1 / 2 . The transit time of plasma 
along their jets leads to a lag of order 10 6 to 10 7 years 
between changes in the central engine being reflected in 
the radio lobe emission. The observed radio luminos¬ 
ity from the cocoon may therefore reflect an earlier time 
with higher accretion rates than the present. We spec¬ 
ulate that these objects predicted by our model may be 
“dying” quasars. 

7.2.2. Thin disks 

The jet production efficiency of a thin disk accretion 
flow AGN (observationally identified with a HERG) is 
similarly derived using Equation 5 of lMeieil (]2001l) . That 
is, 

ctd // 3/2 = 0.005(1 + 1.1.7+0.29/) rh 0 ’ 2 . 

(19) 

The highest efficiency that can be produced by this 
mechanism for the accretion rates of our sample is 
£td// 3 ^ 2 = 0.015. Two of the nine HERGs in our 
sample have lower jet production efficiencies, however, 
seven have efficiencies exceeding this value rising up to 
ejet,HERG // 3//2 = 0.09. These observations may be ex- 
plai ned by a “combined disk” accretion fl ow for rh < 0.01 
(e.g. iNaravan et al.l[l996l : IPu et al.ll2012l noting that the 
high efficiency sources have such accretion rates). These 
accretion flows have a geometrically thin outer disk but a 
puffed up inner region like an ADAF. Black holes fuelled 
by such accretion flows could produce jets moderately 
efficiently due to the ADAF core but still exhibit the 
strong emission lines of HERGs. 

8 . CONCLUSIONS 

We presented a new dynamical model for jet-driven 
radio source evolution. This model differs in three re¬ 
spects from previous models: (1) it includes both the su¬ 
personic and subsonic cocoon expansion phases, (2) uses 
external pressure profiles based on semi-analytic galaxy 
formation models and X-ray observations, and (3) in¬ 
cludes the Rayleigh-Taylor mixing of the cocoon with 
the surrounding material. We applied our model to a 
low-redshift volume-limited sample of radio AGN, and 
estimated the jet powers and active lifetimes of these ob¬ 
jects. Radio sources in massive galaxies were found to 
remain active for longer, spend less time in the quiescent 
phase, and inject more energy into their hosts than their 
less massive counterparts. The majority of AGN energy 
output in the low-redshift universe is provided by sources 
in massive hosts with ~ 10 25 WHz _1 luminosities. The 
relationship between the AGN jet powers and the stel¬ 
lar masses of their host galaxies, Q to t M* 0 ' 1±01 , is 
in agreement with the theoretically expected (lack of) 
dependence for maintenance-mode AGN feedback. The 
host haloes of these AGN are likely to be in or close 
to long-term heating-cooling balance, when considering 
both thermal conduction and the suppression of gas cool¬ 
ing by the AGN. Our dynamical model can also explain 
the observed FR-I/II radio morphology dichotomy. 

High-excitation radio galaxies were found to have ac¬ 
tive lifetimes on average a factor of twelve lower and jet 
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powers a factor of four higher than the low-excitation ra¬ 
dio galaxies. This is consistent with the sporadic nature 
of their merger-driven cold gas accretion. The bimodal¬ 
ity in LERG/HERG accretion rates is consistent with 
the divide between the ADAF and thin disk accretion 
flow mechanisms at approximately m = 0.01. The high- 
excitation radio galaxies in our sample were found to 
be a factor of three times less efficient at producing jets 
than the LERGs. The derived jet production efficiencies 
are broadly consistent with the ranges allowed by radia- 
tively efficient and inefficient flow jet production models. 
The jet production efficiency of the lowest accretion rate 
HERGs exceeds the maximum predicted for a thin disk 
accretion flow, which we suggest provides support for the 
“combined disk” accretion flow model. 

The improved characterisation of the interaction be¬ 
tween AGN and their host galaxies resulting from this 
work should provide valuable input into galaxy formation 
and evolution models. This model will also be a useful 
tool in the interpretation of high redshift survey data, 
especially for Square Kilometre Array (SKA) pathfinder 
surveys. An important development will be extending 
the model to estimate source ages from observed spectral 
indices, since the resolution of the radio cocoons may not 
be possible except for the largest sources (e.g. 80 kpc at 
z ~ 1 for ASKAP EMU). At lower redshifts, future sen¬ 
sitive X-ray surveys of the gas density and temperature 
profiles of the AGN environments (e.g. ATHENA0) will 
be complementary to our work. 
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APPENDIX 

ELLIPSOIDAL TRIGONOMETRY 

The working surface pressure at each angle 9 differs from that along the jet axis due to the non-spherical nature of 
the cocoon. There are two relevant pressures here: the external pressure due to the ambient medium, and ram pressure 
due to the cocoon expansion. The latter requires a conversion from the radial expansion rate of each cocoon element 
to the velocity normal to its surface. Sections of the cocoon expanding supersonically will maintain self-similar growth 
(approximately self-similar when outside strong-shock limit) and can thus be modelled by an ellipsoid with a constant 
axis ratio A. The component of the velocity normal to the surface at angle 6 is related to the expansion rate along 
the jet axis v±(9 = 0) through 


m 


v±(9) 
v±_(9 = 0) 


A 2 sin 2 9 + cos 2 8 
A 4 sin 2 9 + cos 2 9 


(Al) 


where we have defined a dimensionless velocity £. This coefficient can be calculated with a typical value of the cocoon 
axis ratio for powerful FR-II sources. We now need to relate this normal velocity along the jet axis, v±(9 = 0), to 
the rate of cocoon expansion in the 9 direction for the assumed ellipsoidal cocoon. The radial distance to the working 
surface at angle 9 from the central black hole is related to the radius along the jet axis R(9 = 0) by 


R(9) _ R(9) _ 1 

R{9 = 0) V±(9 = 0) yj y)2 gj n 2 Q _|_ cog 2 Q 


(A2) 


where we have defined a dimensionless radius rj. The velocity normal to the surface at each angle 9 is therefore 
related to the radial expansion rate of the relevant cocoon volume element through u_l = CR/rj. 

These definitions for the dimensionless radius and velocity are also useful in describing the initial self-similar growth 
phase of the entire cocoon. Here, these equations allow for analytically tractable solutions for the expansion of each 
volume element and thus the fraction of the jet power required by each element to maintain self-similarity. In the 
strong-shock supersonic limit we find the jet power is proportional to the angular location of the volume element as 


8nk sin 9d9 



/3)R 2 ~PR 3 + 2R 3 ~ p RR 


0=0 


r 1 3 ~P{9)C 2 {8) 


(T c -l)QdA(0), 


3 (T x + 1) 
using Equations E3 El and EO 


(A3) 













